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Abstract: When the semi-positive cosmological constant is dynamical, the naive Eu- 
clidean Einstein action is unbounded from below and the Hartle-Hawking wavefunction of 
the universe is not normalizable. With the inclusion of back-reaction (a crucial point), the 
presence of the metric perturbative modes (as well as matter fields) as a radiation term is 
introduced by quantum fluctuation. They act as the environment (that is, to be integrated 
or traced out), and introduce a correction term that provides a bound to the Euclidean 
action. As a result, the improved wavefunction is normalizable. That is, decoherence plays 
an essential role in the consistency of quantum gravity. In the spontaneous creation of 
the universe, this improved wavefunction allows one to compare the tunneling probabilities 
from absolute nothing (i.e., not even classical spacetime) to various vacua (with different 
large spatial dimensions and different low energy spectra) in the stringy cosmic landscape. 
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1. Introduction 

In the probing of the origin of our universe, a particularly attractive idea is the tunneling 
from absolute nothing (here, nothing means not even classical spacetime) [1], or equiva- 
lently, the Hartle-Hawking (HH) no-boundary wave function of the universe [2]. The basic 
idea is illustrated in Figure 1. It is strongly believed that, in superstring theory, there is a 
vast number of stable and meta-stable vacua, with up to 9 or 10 large spatial dimensions. 
This is the cosmic landscape. If we can reliably calculate the tunneling probability from 
absolute nothing to any point in this vast landscape, one may argue that the origin of our 
universe should be the point in the landscape with the largest tunneling probability. This 
tunneling probability is given by (the absolute square of) the wavefunction of the universe. 
However, the HH wavefunction is not normalizable. Furthermore, it gives an answer that 



- 1 - 



j \ Time 




Nothing 



Figure 1: Tunneling from absolute nothing to a deSitter or Inflationary Universe. 

contradicts observations. It is well-known that the problem lies in the infrared/macroscopic 
limit, so quantum/stringy corrections will not be helpful. In Ref [3], we conjectured that 
a normalizable wavefunction results if decoherence effect is included. One can then apply 
this improved wavefunction to find the preferred vacuum in the cosmic landscape (the one 
emerging from tunneling from nothing with the largest tunneling probability), which then 
evolves to today's universe. We call this proposal the Selection of the Original Universe 
Principle (SOUP); that is, today's universe must lie along the road that starts with the 
original preferred vacuum, arguably the 4-dimensional inflationary universe supported by 
observations. In this paper, we investigate more carefully the decoherence effect on the tun- 
neling probability/wavefunction. As expected, the result supports the basic underpinning 
of SOUP, though the details are somewhat more involved. 

Consider the simple Einstein theory in 4 dimensions with a positive cosmological con- 
stant : 

Se = -J^q J d A x^~g (R - 2A) (1.1) 

The tunneling probability from nothing to a closed deSitter (or inflationary) universe with 
cosmological constant A is given by 

r = |*| 2 ^ e -^ s ^ = -§x (L2) 

where Se is the Euclidean action of the S 4 instanton [5]. Suppose A is dynamical, as in 
a model with 4- form field strengths [6]. The tunneling probability T seems to allow us to 
pick out the universe with the largest probability. However, the Euclidean action Se is 
unbounded from below as A — ► 0, so V — > oo. This has at least 2 obvious problems : 

(1) The wavefunction VP ~ e 3n / 2GA [ s no t normalizable. Furthermore, one can easily 
show that there are other topological instantons with even more negative Euclidean action 
and so larger (actually, infinitely larger) tunneling probability [7] . This implies the presence 
of an inconsistency. 
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(2) Phenomenologically, since r — > oo as A — > 0, it will imply the preference of tunnel- 
ing to a flat universe with zero A, which contradicts the big bang history of our universe. 
Since the size of the universe, the cosmic scale factor a, ~ l/VA, this will imply that 
the tunneling is to a universe of super-macroscopic size, contradicting one's intuition that 
tunneling is a quantum process and so should be microscopic. 

This issue is an outstanding problem since the early 1980s. The inconsistency has 
prevented the proper application of the whole idea. Possible resolution to this problem has 
been suggested (see [4] and references therein): 

• Instead of the usual t — > —it, one may choose instead to rotate time to Euclidean 
time via t — > it [8] . This may work for pure gravity, but the inclusion of a scalar field 
(e.g., as required for inflation) leads to catastrophic consequences, since the scalar 
field theory becomes unbounded from below [9]. 

• One may argue that this problem will be corrected by quantum corrections or string 
theory corrections. However, it is easy to see that this is unlikely to be the case. 
Note that the problem occurs for small A, or large universe, since the cosmic scale 
factor a ~ l/\/A. So this is more like an infrared or macroscopic problem than an 
ultra-violet problem. In fact, one can easily see that the loop correction [10,11] does 
not solve the problem. Also, recent work [12], where the exact HH wavefunction is 
obtained in topological string theory, it seems that A — > is again preferred, just like 
the original HH wavefunction. 

• This last property, namely, that the unboundedness problem becomes acute when 
the deSitter universe becomes large, or macroscopic (actually super-macroscopic), 
naturally suggests that the resolution should lie in decoherence [3] . 

In Ref. [3], we argue that the mini-superspace formulation is inadequate, and we pro- 
pose that the inclusion of decoherence effects due to other modes provides a lower bound 
to Se- We then speculate how the improved wavefunction may be used to select the 
stringy vacuum with the largest tunneling probability from absolute nothing. In this pa- 
per, we shall show that decoherence indeed provides a bound to the Euclidean gravity 
action, though the formula for the tunneling probability may be more involved. Here, we 
shall focus on the pure 4-D Einstein gravity case. Generalization to more dimensions is 
straightforward and will be briefly discussed. 

The basic idea of the approach is widely used in physics. Given any complicated 
problem, we usually follow only a limited set of degrees of freedom, called the system. The 
remaining degrees of freedom, called the environment, are either ignored, or, in a better 
approximation, integrated out. In decoherence, integrating out the environment can cause 
the quantum system to behave like a classical system. In effective field theory in particle 
physics, the massive modes are integrated out to produce higher dimensional operators 
(interaction terms) for the light modes. A famous example is the integrating out of the W 
and Z bosons in the electroweak theory that yields the 4-Fermi weak interactions. Another 
example is the integrating out of the hidden sector in supergravity phenomenology. (Note 
that this has nothing to do with loop corrections.) In the Wilson approach in quantum 
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field theory, where high momentum modes are integrated out, and in the above cases, it 
is clear that the physics may crucially depend on the effects coming from integrating out 
the unobserved modes; that is, they cannot be ignored. Of course, loop corrections can 
also induce dynamics/ interactions that are not present in the tree level, as for example in 
the Coleman- Weinberg model, in light-light scattering, and in the running of couplings. 
In QCD, the running of the coupling emerges from the renormalization group improved 
quantum corrections, not just a normal 1-loop effect. We argue that this is also the situation 
in the study of the wavefunction of the universe. One may consider our result as a back- 
reaction improved quantum correction, not just a normal 1-loop effect. 

The basic idea applied to tunneling is quite simple. It is well-known that the quantum 
tunneling of a particle with mass M, or the system, is suppressed if it interacts with an 
environment. Consider a particle at q = in the potential V(q) as shown in Figure 2. In 
the WKB approximation, its tunneling rate is given by 

T ~ exp(-Sb) 

rlo 

S = / V / 2MV(q)dq (1.3) 
JO 

where So is the Euclidean action of the bounce, i.e., the instanton solution [13]. Note that, 
for V(q) bounded from below, So is bounded from below, as required by consistency. 

In a more realistic situation, the particle interacts with a set of other particles, or 
modes, say x a . However, we are only interested in the quantum status of q, so these 
other modes are integrated out in the path integral, or traced over in the density matrix 
formulation. They provide the environment. Their presence typically introduces a frictional 
force to the evolution of q. It was shown by Sethna [14] and by Caldeira and Leggett [15], 
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that the bounce So increases to S ~ So (1 + fj) where f] > is proportional to the coefficient 
of friction (see Appendix A for details). That is, the interaction with the environment 
suppresses the tunneling rate. (This suppression takes place even if no friction is generated.) 
One may understand this result in a number of (equivalent) ways : 

• As a quantum system, 



That is, the increase in S is due entirely to the longer path length in the many 
dimensional (q, x a ) space. 

• The interaction of q with the x a interferes with its attempt in tunneling. One may 
view the interaction with x a as attempts to observe q. Repeated measurements of q 
or repeated attempts of measuring q suppresses the tunneling rate. This is analogous 
to the Zeno or Watch Pot effect. 

• The interaction of q with the environment x a diminishes the quantum coherence. As 
a consequence, the system behaves more like a classical system than like a quantum 
system. Since tunneling is a pure quantum phenomenon, it should be suppressed as 
the system becomes more macroscopic/classical. We shall refer to this as decoherence, 
by which we mean the process where a quantum system behaves more classically (i.e., 
less quantum) via its interactions with the environment [16]. 

Here, we study this effect in quantum gravity, in the tunneling from nothing scenario. 
In this case, the cosmic scale factor a plays the role of the system, while the metric fluc- 
tuations around a (and any matter field modes) play the role of x a , i.e., the environment. 
Figure 3 illustrates this situation. Since we measure only a, the metric fluctuations are 
integrated out in the path integral (or traced out in the density matrix). As expected, we 
shall show that their presence suppresses the tunneling from nothing to a deSitter universe. 
As expected, the decoherence effect is negligible for large cosmological constant (small size 
universe), but becomes increasingly important as A — ► 0, as it suppresses the tunneling 
probability when a = l/yA/3 becomes macroscopic. In contrast to a normal quantum 
system, the inclusion of the environment is of fundamental importance in quantum gravity, 
since the corrected Euclidean action is now bounded from below. 

The coupling of the metric fluctuations and scalar fields to a(r) are more complicated 
than that of the above quantum system. Each metric perturbative mode behaves like a 
simple harmonic oscillator but with time-dependent (or a-dependent) mass and frequency. 
However, the real subtlety of the calculation comes in another way. If we treat the metric 
fluctuation modes as pure perturbations, we shall get nothing except the loop correction, 
a known result in Euclidean gravity. This is not hard to see. The Euclidean action in 
minisuperspace (that is, keeping only the cosmic scale factor) is given by 




(1.4) 




(1.5) 



- 5 - 




Figure 3: Typical tunneling involves metric and matter excitations as well. Loop effects may also 
be included, though they are not important for our discussions. 



where a rescaling has rendered r, a(r), the Hubble constant H = \/^/~\ and A = 2GA/97T 
dimensionless. For S 4 , Ha(r) = sin(ifr), where south pole (north pole) corresponds to 
Ht = (ir). This path gives the Euclidean action (|1.2| ). For a fluctuation mode / that 
satisfies the classical equation, the Euclidean classical action can be written as a surface 
term Se ~ a 3 / f\^ H = 0, since a(r) = at the two poles. As a result, no decoherence 
term is generated. (Including an additional boundary term makes no difference.) However, 
we find that the back-reaction is crucial for getting the correct answer. Instead of using 
the unperturbed a(r) given above for the S geometry, we leave it arbitrary during the 
tracing out of the fluctuation modes. In the path-integral formalism, one starts with the 
path integral that includes the scale factor a as well as the perturbations around a. The 
following, for example, shows the inclusion of the metric tensor perturbations t n 



Z= f D[a] TT f " D[t n ]e- S * exp ( -~ f dra 3 V[t n 2 + ^n(a) 2 ^] I (1.6) 
Tracing out the perturbations, we have 

Tr[Z] = II / < / " tf n) Z (1-7) 

This results in a new term in the modified action 

S E ,dC -\J dT (-°a 2 - « + Aa 3 + (1.8) 

where the last term, coming from integrating out the perturbative modes, behaves like 
ordinary radiation, v is a constant that measures the number of perturbative modes. We 
then solve for a(r) and obtain the corrected Euclidean action Se,cLC m the saddle-point 
approximation. It turns out that the v term modifies the shape of the 5" 4 instanton to 
barrel-like, as shown in Figure 4. Since the effect is perturbative in nature, we expect the 
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Figure 4: Modification of the S instanton due to the presence of decoherence, i.e., the inclusion 
of the environment. 

barrel to have the same topology as S* 4 . We see that, due to the back-reaction, a(r) does 
not vanish at the two ends. However, the contribution of the end plates of the barrel to 
Se,cIC happens to be zero. 

After tracing over the metric fluctuations in this way, the rate of tunneling from nothing 
(i.e., no classical spacetime) to a deSitter universe (much like the inflationary universe) is 
now given by T ~ exp(F) = exp(— Se jC ic), where 



where D is the decoherence term and c depends on the cut-off. In string theory, that cut-off 
is naturally provided by the string scale. Note that Se,cLC is now bounded from below. See 
Figure 5. Here, c in string theory also depends on the string spectrum, so c should be 
calculated for each vacuum. We find that 



where l s is the string scale (a 1 = (l s /2ir) 2 ) and ridof is the number of light degrees of 
freedom included in the environment. For the pure gravity case, we have n^of = 2 for the 
two tensor modes. For small A, below a critical value, the tunneling is actually totally 
suppressed. 

This is quite understandable. A precise S 4 spherical geometry leads to the pure de 
Sitter space, which by definition excludes any radiation. If we do not allow back-reaction, 
then the system a(r) actually cannot feel the presence of the environemnt. Allowing 
back-reaction, the quantum fluctuation during the spontaneous creation of the universe 
generates some radiation (even though no radiation is introduced classically). They act as 
the environment. The presence of this environment modifies the geometry in a way so that 
tunneling is to a universe with both a cosmological constant and some radiation, with a 
suppressed tunneling rate. The presence of the radiation generates the decoherence term. 
Note that, in integrating out the perturbative modes, the case with zero amplitude for the 
perturbative modes is also included. As we see, the amount of radiation is proportional to 
1/A 2 , so there is more radiation in a larger universe. 




(1.9) 




(1.10) 
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Figure 5: Semc and Se.o as functions of the cosmological constant A. We see that, with the 
inclusion of the decoherence effect, the Euclidean action for the S 4 instanton is now bounded from 
below. For large A, the semi-classical approximation breaks down. For small A, tunneling is actually 
completely suppressed. 

What we consider fundamental is that the decoherence effect actually provides the 
Euclidean action Se of pure gravity with a lower bound. Since the metric fluctuation 
contribution to Se cannot be turned off, they must be included in the evaluation of the 
tunneling rate. For usual quantum system, So is bounded from below. So one may view 
this friction/environment /decoherence effect as a correction, albeit it may be very big. In 
quantum gravity, this effect resolves the unboundedness problem. That the quantum fluc- 
tuation provides a natural source to cure the boundedness problem implies that quantum 
gravity is actually self-consistent in the macroscopic regime. 

This effect renders the system more macroscopic and so less quantum. That is, as a 
becomes large, its interaction with the metric fluctuations and matter fields should sup- 
press the tunneling rate. Indeed, we show that the inclusion of the metric fluctuation 
decreases the tunneling rate, as expected. Note that the decoherence term is not the usual 
perturbative quantum correction. For large A, where quantum correction is expected to 
be large, the decoherence term actually becomes negligible. One can include the quantum 
corrections; however, they do not change the qualitative behavior for moderate values of 
A, i.e., GA « 1. 

Let us gain some idea of the magnitude of F. Following Eq.( |1.9[) , we find the value of 
A with maximum tunneling probability is 



For the string scale a few orders of magnitude smaller than Mpi, we can easily have 
F — lO 1 ^ or larger. On the other hand, the critical value of A is A c = A max /2. For A < A c , 
the barrel-shaped instanton is destroyed. At A c , F(A C ) = 0, so the tunneling probability 
at A close to A c is already negligibly small compared to that at A max . 
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Having resolved the outstanding problem mentioned above, one may then apply this 
consistent tunneling approach to select a preferred vacuum in the cosmic landscape in 
string theory It is straightforward to generalize Eq,([L.9|) to arbitrary spatial dimensions. 
In particular, for ten dimensional spacetime, we have 

- F = S E4C * S e ,iq + c ( ^j? ) (1.12) 



710 

where Se,io is the 10-D Euclidean action determined in mini-superspace and V%o is the 10- 
dimensional volume of the instanton. Note that (see Ref [3]) Se.io reduces to —3tt/GA for a 
vacuum where the extra 6 dimensions are compactified. For each vacuum, c depends on the 
spectrum. It may also depend on the compactification and dilaton moduli. Knowing these 
properties of each vacuum allows one to calculate its tunneling probability from nothing. 
One may estimate the size of c by comparing to the 4-D case: 

n ^ ] — (1.13) 



with 8irG = 1/Mp t and g s is the string coupling. Since the string scale is expected to be a 
few orders of magnitude below the Planck scale, we expect c to be a small number. In [3], 
we consider the suppression of tunnelling in the context of the spontaneous creation of the 
universe. The suppression happens due to the effect of both the gravitational and matter 
perturbations on the bounce solution to the Euclidean Einstein equations. The general 
idea, loosely speaking, is to separate the "universe" into a "system" (the pure gravita- 
tional bounce) and the "environment" (the perturbations). One is interested in measuring 
properties of the system and in order to do so one simply traces out the environmental 
degrees of freedom. The effect of the environment can be significant and such effects have 
been studied in various systems. In the previous work [3] we estimated this effect on the 
tunnelling probability of the universe by considering the unperturbed deSitter space as the 
"system" and the "environment" consists of the metric perturbations and matter fields. 
We see that the qualitative features are robust, though the details are somewhat more 
involved. (There we used Vg instead of V\o for Eq.( 1.12]) ). For large A, the semi-classical 



approximation breaks down. For small A, additional decoherence effect may be important. 
(Here we have calculated only the leading order.) Fortunately, the range of A where the 
tunneling probability is largest seem to lie in the region where the approximation is most 
reliable; and we are interested in vacua with large tunneling probabilities. Since the tun- 
neling probability drops off rapidly as A — > and A — > oo, their precise values are not as 
important to us. With some luck, we may use the above formula to locate the preferred 
set of vacua in the cosmic string landscape. Note that c depends on the spectrum at each 
point in the landscape. 

As one can see in Figure 5, intermediate values of A, much like the inflationary universe 
that describes the history of our universe [17] seems to be preferred. As in Ref [3], we find 
phenomenologically that 10-dimensional deSitter-like vacua (with F ~ 10 9 for instantons 
S 10 , S 5 x S 5 , S 4 x S 3 x 5 3 ) are not preferred while super symmetric vacua in any dimension 
have essentially zero tunneling probability. Also, tunneling to vacua very much like our 
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today's universe (with a very small dark energy) seems to be severely suppressed. Tunneling 
to a universe with quantum foam seems not preferred either. Among the known vacua, the 
preferred ones are the 4-D brane inflation [18] as realized in a realistic string model [19], very 
much like the inflationary universe that our universe has gone through (with F ~ 10 16 ). 
Although the details of the decoherence term obtained here is somewhat different from 
that used in Ref [3], we see that the qualitative features summarized there remain true. 

Decoherence and related issues in quantum/Euclidean gravity have been studied earlier 
[20-24], where they are mostly concerned with the evolution of the inflationary universe. 
Here, we study decoherence in the quantum tunneling in gravity. To be self-contained, we 
shall review some of the relevant formalisms developed there. 

In Sec. 2, we review the tunneling to closed deSitter space and perturbative modes 
around the cosmic scale factor a(r). In Sec. 3, we evaluate the contribution due to 
these modes to the effective Euclidean action. Here we find that no decoherence term 
is generated in the S 4 background. To see why back-reaction is important, the effect is 
calculated in the background of a "squashed" S 4 geometry. The reader may choose to 
skip this section. In Sec. 4, we study the modified bounce including back reaction. We 
derive the effective Euclidean action after integrating out the metric perturbation with 
an arbitrary background a(r). In Sec. 5, We find the modified bounce solution from 
the effective Euclidean action. This allows us to obtain the improved wavefunction and 
tunneling probability from nothing. In Sec. 6, we show that the new term in the improved 
Euclidean action behaves like ordinary radiation. We see also how the new term appears 
in the Wheeler-DeWitt equation and how it influences the tunneling amplitude. In Sec. 7, 
we discuss the implication of the main result and its connection to the proposal of Ref [3]. 
Sec. 8 contains some discussions and Sec. 9 contains a summary and further remarks. 
Some of the details are relegated to appendices. For the sake of completeness, a number of 
known results are reviewed extensively. 

2. Setup 

2.1 Notations 

We shall render various quantities dimensionless for the ease of calculation. The conventions 
we follow are those of [2]. The Euclidean action is defined as 




(2.1) 



The Euclidean metric is given by 



ds 2 = 



a' 



2 (dr 2 + a(r) 



2 dn 2 ) 



(2.2) 



where is a 



-2 



2G/37T. With this metric ansatz, the action becomes 




(2.3) 



where A = cr 2 A/3 = 2GA/9vr. So A, r, and a are all dimensionless. 
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2.2 The deSitter Space 

In this section we give a summary of the result of [20]. Consider a compact three-surface 
S 3 which divides the four-manifold M into two parts. One can introduce the coordinates 
x % (i = 1,2,3) and a coordinate t such that S 3 is the surface at t = 0. The metric takes 
the form 

ds 2 = -(N 2 - NiN^dt 2 + 2N i dx i dt + h ij dx i dx j . (2.4) 

where N and TVj are the lapse function and the shift vector, respectively. The action is 
given by 

S = j " (L g + L m ) d A x (2.5) 

where 

L g = (G^KijKu + VhR) (2.6) 

is the Ricci scalar of t 

given by 



where R is the Ricci scalar of the three-surface and is the second fundamental form 



In the above expression "|" denotes the covariant derivative. G^ kl is called the metric of 
the "superspace" and is given by 

G ijkl = 1 ^ h ik h ji + h u h jk _ (2.8) 

In the case of a massive scalar field, the matter Lagrangian L m is given by 

1 / /<9<I>\ 2 N { d$ d$ 

- (h« - - m 2 * 2 

V J dxi dxi 

In the Hamiltonian treatment of general relativity one treats hij and $ as the canonical 
coordinates. The canonically conjugate momenta are 

dL„ VhM? 



7T 



ij — 9 _ v "> lvJ -p 

dhij " 16vr 
dL m 

7T<J> = — 

9$ 



(K ij - h ij K) (2.10) 



V dx l J 



The Hamiltonian is 



H = j d 3 x(n ij h ij +Tr< s> <&-Lg-L m ) (2.11) 
= j d 3 x (NHo + NiW) 
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where 



M 2 

H = l6i:Mp 2 G ijH 7r^TT kl - -^^hR+ (2.12) 

l07T 



ax* ax-? 



^ = -2^ + ^^ (2.13) 



and 



Gijki = -h x l 2 (h ik hji + huhj k - hijh k i) (2-14) 

The quantities TV and TVj are regarded as Lagrange multipliers. Thus the solution obeys 
the momentum constraint 



and the Hamiltonian constraint 



2.3 The Perturbations 



W = (2.15) 



H = 0. (2.16) 



Now we study the perturbations around the deSitter spacetime. The perturbed deSitter 



has a three-metric of the form 



hij = a z (SUj + dj) (2.17) 

where fiy is the metric on the unit three-sphere and is a perturbation on this metric 
and can be expanded in harmonics: 



n.l.m 



^ c nlm{Sij)?m + ^ c nlm(^ij)?m + 2 ^ m (G°)JJ n + 2t^ m (G^)™ m (2.18) 

where the coefficients a n i m , b n i m , c° nlm , c e nlm , t^ nlm and t e nlm are functions of time but not 
the three space coordinates. The Q(x l ) are the scalar harmonics on the three-sphere. The 
Pij(x l ) are given by 

Pij = tf—[Q\ij + 3 n «Q ( 2 - 19 ) 
where we have suppressed the n,l,m indices. The Sij are given by 

Sij = + Sj\i (2.20) 
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where Si are the transverse vector harmonics. The Gij are the transverse traceless tensor 



harmonics. Further details can be found in [20,25]. The lapse, shift, and the scalar field 
&(x l ,t) can be expanded in terms of harmonics as 



N = Ni 



1 + ^2 9nlmQ 



n.l.m 



Ni = a(t) Y, 



n,l,m 



n,Lm 



where Pi = ^TZ\Q\i- The perturbed action is now given by 

S = S M,N )+YSn 



(2.21) 



(2.22) 



where we have denoted the labels n, I, m, o, and e by the single label n. Sq is the action 
of the unperturbed deSitter, 



If o / a 2 1 rh* 



(2.23) 



S n is quadratic in perturbations and is given by 



S n = J dt{L n g + U r 



(2.24) 



where L™ and are the nth mode gravitational and the matter Lagrangians, respectively. 
As we are interested only in the gravitational tensor and the scalar field perturbations in 
this paper, we only display those terms here. In the absence of sources, the other modes 
can be gauged away. For further details and the intricacies of all perturbations we refer 
the reader to [20]. 

The tensor perturbations t n have the Euclidean action 



&n = 2 / tnDtn + boundary term 



1 



where 



D 



d_ 
"ah 



dr 



+ a(n 2 - 1) 



(2.25) 



(2.26) 



where the background satisfied the classical equation of motion. The action is extremized 
when t n satisfies the equation Dt n = 0. Setting Nq = 1, this is just the equation 



d_ 
~dt 



[a 3 t n + (n 2 - l)at n ] = 



(2.27) 
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For t n that satisfies the equation of motion, the action is just the boundary term 

S^ = \a"{t n i n + ^ (2.28) 
The path integral over t n will be 

J d[t n ] exp(-Sf ( d )) = (detD)- 1/2 exp(-S 



E(cl) } 



(2.29) 



Scalar fields can be treated in a similar fashion. The scalar field perturbation La- 
grangian is given by 



1 



N a 3 



1 



fn m f n „o fn 



(2.30) 



where we work in an appropriate gauge choice. Setting iVo = 1, the equation of motion for 
fn is 



d_ 
~dt 



f n ) + [m 2 a 3 + (n 2 - l)o] f n = 



(2.31) 



One can evaluate the scalar field path integral and the results are similar to that of the 
tensor perturbative modes. So in this paper we deal only with the gravitational tensor 
modes. 



3. Perturbative Correction to the Bounce : No Back Reaction 

Before doing the calculation that includes the backreaction due to the perturbative modes, 
we first perform a simple calculation to illustrate the key issues we are facing. The result 
of this section is meant to be a warm up and indicates the possibility of a correction to the 
wavefunction. We consider a 5 4 solution with both its polar regions flattened by a small 
parameter 5 and show that this squashed geometry allows for extra perturbative modes 
that can have significant effect on the calculation of the wavefunction. When 5 = 0, these 
extra perturbative modes vanish and what remains is the well known one-loop correction to 
the S 4 . For a proper treatment the reader may go directly to the next section and onward. 

The Euclidean equation of motion for the nth tensor perturbation mode on a S 4 
background follows from Eq.( 2.27| ) 



do ( (n 2 — 1) \ , 
t n + 3—t n -(± —L)t n = (3.1) 

where dot denotes differentiation with respect to the r variable and a (r) = A -1 / 2 sin(y/~\T). 
This can be converted to a more familiar form by the substitution F n = \f\a {r)t n . In 
terms of F n and x = cos (At) the equation of motion becomes 

9x d 2 F„ dF n ( n 2 \ , , 



dx 2 dx \ (1 — x 
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which is just an associate Legendre equation of degree one and order n. This has two 
linearly independent solutions, P 1 ~ n (x) and Q™(x). However, because a Q = at x = ±1, 
we have 

x=-l 

= (3.3) 

x=l 

We shall see in the following sections that once the perturbative modes are prop- 
erly accounted for, S 4 changes to a "barrel". Anticipating this result we consider the 
S instanton that is flattened slightly at the two poles. This "squashed" S* 4 is given by 
a Q (r) = ^ sin(Ar), with the restricted range 1 — 5 > x > —1 + 5. There is a good reason 
for considering the squashed geometry. The perturbative modes can be strong enough to 
change the geometry of S 4 (and as we shall see, they will). So fixing the geometry to S 4 
before doing the perturbative analysis is too restrictive. The perturbative modes on S 
have to vanish at the two poles if they are to respect the background geometry. Their effect 
has been studied in [11,32] and apart from contributing to one-loop effect they do not lead 
to any tunneling suppression. However, to get the decoherence effect one must include 
modes that have nonvanishing values at the two poles. And we shall see that the squashed 
geometry allows for modes that can potentially lead to decoherence. This expectation will 
be confirmed in the later sections. 



SI 



i \3j. dt n 

CLT 



As explained in the introduction (Eq. (|1.6 ,1.7)), to trace out a given perturbative mode 
(say, t n ) one must perform a path integral over that mode with the initial and final ampli- 
tudes the same (say, t l n = tn), and then integrate over all possible values of t n . This is just 
the trace operation in the path integral formalism. To do this we first find the action for 
t n (details can be found in the appendix) 



1 



I \3j. 

dr 



2 

5(2 - 5) 



-|X=(-l+<5) 
x=(l-<5) 



(3.4) 



(1-5) + 5(2 -5) 



d x Q^(l~5) 
Q^l-5) 



(t 



i \2 



Next we must find the prefactor to the path integral for t n . This can be found as explained 
in the Appendix B. The prefactor is given by 



2 J 



/4vr 2 n 
Integrating over all initial states t n (tracing out the nth mode) gives 

dt n K (4 x = -1 + 5; 4 x = 1 - 5) = a(n)f(5)2- n 



(3.5) 



(3.6) 



where a(n) contains factors of n and f(5) is a function of 5, with f(0) = 0. K (t l n , — 1 + 6; t n , 1 — 6) 
is the path integral over t n mode with the boundary conditions t n (x = — 1 + 5) = t n (x = 
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1 — 5) =t l n . Taking care of all the modes by tracing over all of them (within the lower and 
upper cut-offs), we get 

[] / <tf(C-l + 5;4l-<S) =l\a(n)f(5)2- n (3.7) 

n J n 

= A(n,5)e~ N ^ 21 ^ 

where N counts the modes. As explained later, 

where v = g^-^r- The wavefunction then corresponds to the path integral over the scale 
factor a and the perturbative modes t n with the full action as given in Eq.|E22D. We trace 
over the perturbative modes and get 

$ ~ e (jx- D ) (3.9) 

where D is the decoherence term leading to tunneling suppression 

D ~ ^ (3.10) 

so the above result leads to an order 0(1/A 2 ) suppression to the Hartle-Hawking wavefunc- 
tion. This naive (naive because, as it turns out, the backreaction is important and also 
leads to an 0(1/A 2 ) contribution) expectation is indeed vindicated in our calculation of the 
modified bounce. The more careful calculation changes the coefficient hi(2), but maintains 
the inverse square dependence on A. 

4. The Modified Bounce : Including Backreaction 

The bounce solution S 4 is a solution to the Euclidean Einstein equation which is obtained 
as the Euler-Lagrange equation from the Euclidean action Sq [a]. 

a d 2 1 , , , 

-2 + =3A (4.1) 

a o z a z 

The solution is given by a(r) = y^' cos (v / At). This is the bounce in the absence of any 
perturbation. Including the perturbations, treating them as the environment, and tracing 
them out (as explained in Eq.(1.6, L7|) ) will lead to a modified equation of motion instead 



of Eq. (|4.l|) . In this section we derive this modified bounce equation. 

Let us consider the effect of the metric perturbations. To find the modified bounce 
equation we have to carry out the path integral over the perturbation modes t n is Eq. 



( 2.29 ). Let us write down the path integral for a single tensor perturbation mode 

D[t n (r)} exp ( -i [ T/2 dra 3 [t n 2 + ^^i 2 ] ] (4.2) 



/ 
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This is a path integral for an oscillator with a varying mass as well as frequency. We can 
simplify this to a path integral of an oscillator with constant mass and variable frequency 
(given by uj n {a(u)) = \J (n 2 — l)a(u) 2 ) using a new variable u 

The Euclidean action now becomes 

I r u f 

Se = ^ du{t'Z + {n 2 -l)aH 2 n ) (4.4) 

To keep notation uncluttered, we drop the subscript n for the time being. Let t = t c i + 1, 
where t c i is a solution to the classical equation of motion for the above action 

t" d - io n (a) 2 t d = (4.5) 

with t c [(ui) = t(ui) and t c i(uf) = t(uj). Here t denotes fluctuations about the classical 
solution with t(ui) = i(uf) = 0. That is, (tt')\u( = (t c ii')\u{ + (*d*d')|uf ■ Apriori, the 
contributions to the path integral will come from t c / and t. Here, the fluctuations t will lead 
to a prefactor. We shall keep track of this prefactor as it will have important contribution 
to the modified action. Substituting t = t c [ + i in Eq.( [4.4[ ), we obtain 



Se = \ du {-it" + u 2 n P) + {t cl t cl ')\ u ui (4.6) 



The second term in the above equation is simply 5^(t c /). The path integral in Eq. ([4.2|) is, 
therefore, given by 



e -s E (t cl ) p[i]=e ~s E (t cl ) f D [i( u )]ex V (-l r du(-ti" + ulP] 

Jo V 2 J u 



(4.7) 



Notice that the prefactor P[t\ is independent of t(m) or t(uj). So the reduced path integral 
Eq.(|i~7^) becomes 

Tr[Z] = I D[a\l[P[t n ] [ dt n ! dt&fa - ti)e~ S ^) (4.8) 

We shall first evaluate the integral over t l n = td. In Appenxix C we evaluate the prefactor 
P[t n ]. It is straightforward to include a scalar field. 

Although finding a general solution to Eq.( f4.5| ) is in general impossible, we note that 
it has the same form as the Schrodinger equation for a particle in a potential uo n {a{u)). 
So for a slowly varying potential we can use the WKB method for finding the solution to 
Eq.([4.5|). By slowly varying one means that the following condition is satisfied 
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This is satisfed by the higher modes n >> 1 and we shall see that these modes are the 
ones that contribute to the suppression of quantum tunneling. (One may be concerned 
that uj n —* as a — > 0. As we shall see, in contrast to the unperturbed bounce, a actually 
stays finite in the modified bounce. In this sense, the WKB approximation is reasonable.) 
The two independent WKB solutions to Eq.([4,5|) are 

t± = -^=exp(±J du'oj n (u')^ (4.10) 

The general solution will be a linear combination of the independent solutions 

t cl (u) = At+(u) + Bt d (u) (4.11) 

The values of A and B will depend on the boundary conditions. If Eq.(4.11) is to satisfy 
the following boundary conditions 

t cl {ui) = t l t cl {u f )=tf (4.12) 
then we have the following values of A and B 

a _ t f ^u n {u f ) - t % y/u} n (uj) exp(-D n ) 

(exp(A0 - exp(-AO) ' 

t f s/u^Juf) - t % ^LO n (ui) exp(£>„) 
a = — 



(exp(D n ) - exp(-D n )) 



where D n is given by 

D n = I du io n {a(u)) = I dr " W 7 V ' 3 " (4.14) 

Next, substituting the solution Eq.( 4.11| ), using Eq.( |B.2D , in Eq.(4.6), we get 

S E (t d ) = \ (t cl {u f ) d ^l - t d { Ui )^-) (4.15) 

= 2(exp(£>„)-cxp(-g„)) Wf^n{Uf) + (j*) 2 U n (Ui)) (exp(D„) + exp(-D n )) 

-4tH f y/u n (ui)uj n (uf)\ 

This gives the classical contribution Eq.( |4.7D . Note that in order to perform the trace over 
this perturbation mode, we will be setting t l = t* , and u> n (ui) = u> n (uj). The path integral 
over t leads to the prefactor. Its evaluation for a time-dependent frequency is explained in 
Appendix C (Eq. fl02| )). (Equivalently, one may use the Pauli-van Vleck-Morette formula 
studied by Barvinsky [26]). The prefactor (path integral over t) gives 



D[i(u)} exp (\ r du (-& + ulP)) = ( f^ i) ( =~r (4-16) 
\ 2 Ju, J V(exp(£y) - exp(-Ai)) 

Now the trace over the perturbative modes can be performed by setting t % = v and 
integrating over the amplitudes t l 

f dtf f dPSfi - tf) D[t(u)\ exp (S E ) (4.17) 



1 



/ dtf exp(-S E [td\) 



yj (exp(D„)-cxp(-D n )) 
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where we used f dt$ exp (— SE[t c i\) — l/-v/w n («/). Also, note that D n » 1 for large 
values of n, and (exp(Z? ra ) — exp(— D n )) -1 ' 2 ~ exp(— D n /2) . Finally, one has to perform 
this tracing-over over all n modes. This leads to the following contribution 

n n twt i m " exp( - s ( 4 - i§ ) 



From the expansion in Eq. fl2.lq ) it is clear that one has to count the indices n,l,m (cor- 



responding to the spherical harmonics on S s ). For a given n, I can take n — 2 values 
(Z = 2, 3, n — 1), and for a given I, m can take 2Z + 1 values {m = —I, —I + 1, I — 1, Z). 
For a given n this introduces a degeneracy of f(n) ~ 2n 2 . The factor of 2 comes due to 
the two tensor modes t^ lm and t^ Zm in the expansion in Eq.( 2.18| ). To proceed further, we 
need to define the cut-offs. A natural short wavelength cut-off is the string scale, l s . The 
long wavelength cut-off is the inverse Hubble length. This gives a cut-off for n 

— -"-(T^-ih/s (419) 

where we have used Eq.(4.14). Note that uj n (a) 2 = (n 2 — l)a 4 . For higher modes, therefore, 
u n ~ na 2 . Thus, 

1 rT/2 dr ^ 

D = - — > n (4.21) 

2 7-T/2 a V 

' n,l,m 

N 4 f T ' 2 dr 1 /- T / 2 , !/ 

ar 

f_T/2 « z V-T/2 



(X 2 /_7 n /9 A 2 a 



We can also easily generalize to arbitrary number of degrees of freedom, with 

v _ n do fN 4 n do f 



(4.22) 



A 2 4 4(tfZ s ) 4 

For pure gravity, we have n do j = 2. The modifed Euclidean action is, therefore, given by 

S E ,dC = S E ,o + D[a] = \f dr (-ah 2 - a + Aa 3 + (4.23) 

To evaluate Se dC via the saddle-point approximation, we have to find the classical path, 
or the bounce solution, of this effective action. 

5. The Bounce Solution 



The equation of motion for the Euclidean action, Eq.( [4.23 ), is 
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Figure 6: The Modified Bounce : It is a squashed version of S 4 resulting in a "barrel". 



This is just the Euclidean version of the Einstein equation with both a cosmological constant 
and radiation. As is well known, this equation allows for a variety of solutions (Ml, M2, 
Al, A2, E, 01) [27]. One has to be careful in deciding which one of these is the correct 
bounce solution. The bounce solution is the real solution to the classical Euclidean equation 
of motion. The M2 solution satisfies the bounce criteria and is given by (in the Euclidean 
form) 




This solution is the modified bounce and has some interesting features. 

• It is a deformation of the S 4 instanton from spherical to barrel-shaped, with — 7r/2\/A < 
r < 7r/2\/A. It reduces to the usual S A instanton when v = 0, as expected 

a(r) ->■ -^cos(VAt) (5.3) 
v A 

Here v is the deformation parameter. 

• For large values of v the bounce is destroyed. Presumably, there is no more quantum 
tunneling because of excessive decoherence. This critical value of A is given by A c = 
\v. For tunneling A > A c . 
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Figure 7: The solid curve is F — —SE.dc, using the exact value of Q\ The dotted curve shows F 
using 3 ~ 8/3. The light flat curve is the integral 9, which asymptotes to 8/3 for large A. The 
plots are given with Mpi/M s = 10 3 . The dotted and the solid curves differ only slightly, justifying 
the approximation made. 



The modified Euclidean action can now be calculated. One has to be careful though to 
include the 0(u/X 2 ) contribution from the —act 2 — a + Aa 3 part. This contribution will add 

V 

A 2 a(r) 



to the 0(v/\ 2 ) contribution from the ^ —> term and give the total 0(u/X 2 ) modification. 



We state the final result 



where 9 is an integral given by 



_f dr (^l) (5.5) 

27K 



o(r) 



where a(r) is the bounce solution given by Eq^p^). The exact value of 9 involves an 



elliptic integral, but we can make a fairly accurate estimate and the result is 

g 

9 ~ - (5.6) 

For tunneling from nothing, we should consider the barrel to be a deformed 5 4 with 
the same topology. That is, we should include the contribution of end plates of the barrel 
in the evaluation of SE,dc i n Eq.(|5Ti|). At the end plates, 

1 / 




2A y 

However, their contribution to SE,dc is zero because the weighing factor sin m 
Eq.(^l) vanishes at the end plates. 
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As expected, for v = we get the usual Hartle-Hawking wavefunction. Now, including 
the environment, the tunneling probability is given by 

p~ e UA (5.7) 

For A — > 0, higher order decoherence effects should be included. These effects can come 
from a more careful treatment of the leading term we have obtained, or from higher order 
interaction of the modes with a and among themselves. Quantum effects may become 
relevant. This is especially so at large A. Note that SE,dc peaks around A = 8v. So for 
regions with non-zero tunneling probability, A > A c = Av. 



6. The Modified Hamilton Constraint and the Wheeler-DeWitt Equation 



The modified action, as given in Eq.(4.23), leads to a modified Hamiltonian constraint. 



The modified Lorentzian action is given by 



S = - j dr (-aa 2 + a - Aa 3 - -^-) (6.1) 
2 / \ X z aJ 



The modified Hamiltonian constraint is 

H = j-(-n 2 a -a 2 + \a A + ^) =0 (6.2) 

where II a = —ad is the conjugate momentum. Using the Hamiltonian constraint and the 
equation of motion following from the modified Lorentzian action, we have 

,d.o 1 v , . 

a v 
a = A^4 

Since a/a = — YliPi + 3pi)/2 = — Yl PiO- + 3^j)/2, this implies that the new v term has 
equation of state uj = 1/3, precisely that of ordinary radiation. 

One can write down a differential equation describing the evolution of the wavefunction 
- the Wheeler-DeWitt (WdW) equation -by imposing the condition as an operator equation. 
This is quantum gravitational analog of the Schrodinger equation in quantum mechanics. 

-M = (^- + U(a) ) J ^ = 

where n a = —id /da and U(a) is the gravitational potential given by 

U(a) = U (a) - ^ = a 2 - Aa 4 - ^ (6.4) 

Recall the case without the v term. In the classically forbidden region, i.e., the under- 
barrier region A^ 1 / 2 > a > 0, the WKB solutions for the tunneling amplitude from a = 
to a = A -1 / 2 are 



; ±/ 1/v ^|n(a')|da' (g_ 5 ) 
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Note the Hartle-Hawking no boundary prescription requires us to take the positive sign in 
the exponent in Eq.( |6.5| ), that is, \&_|_. This yields the HH wavefunction. Including the v 
term is like solving the wave equation with a positive (instead of zero) energy eigenvalue 
E = vj\ 2 (see Figure (§)), 

/rf 2 \ 

*(o) = EV(a) 

The presence of E > decreases the tunneling amplitude Vl/+. To see the origin of E, 
consider the presence of a generic field Xi so the WdW equation is crudely given by 

{~\- U o(°) + ^ + "(X, aA ¥(a, x) = (6.6) 

Note that the kinetic terms in Eq.( |6.6| ) have opposite signs. Let 

*(o,x) = ^ n ^n{a)u n (x) (6.7) 
then the above equation can be separated, 

1 X + w(x, a) 2 Wn(x) = £n(a)n„(x) (6.8) 



2 

1 d 2 Vv, 



+ Uo{a)ij} n = e n (a)if), 



2 da 2 

Taking the ground state, we obtain the zero point energy eo, with ^f(a) = ipo(a)- Including 
the zero point energy of all x fields then yields E, which turns out to be independent of a. 

If instead, we take \P_ in Eq. (|6.5| ), as suggested by Linde and Vilenkin [8] and which 
is more familiar in quantum mechanics, we see that the inclusion of the v term enhances 
tunneling. This means interaction with the environment enhances the tunneling amplitude, 
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which is counter-intuitive. Furthermore, the highest tunneling probability will go for large 
A, when the semi-classical approximation used here breaks down. We believe the Hartle- 
Hawking no boundary prescription is the correct one for the spontaneous creation of the 
universe. 

As we see clearly now, this v term plays the role of radiation. The deformed S* 4 has, 
in a sense, allow the spontaneous creation of a universe with some radiation in it. In 
pure gravity, this is simply gravitational radiation. Due to the presence of such radiation, 
the cosmic scale factor a(r) does not vanish any more at the poles of S 4 . In fact, the 
decoherence has led to a flattening of the two poles. 

We note that, in theories with dynamical A, there will be momentum terms associated 
with A in the Wheeler-DeWitt equation. For example, if A is associated with the potential 
energy of some scalar field (ft, then there will be momentum term in the Wheeler-DeWitt 
equation. However, the inclusion of such a momentum term for A is not expected to change 
the boundedness of the gravitational potential. 

For large cosmological constant A, the radiation is suppressed. In this case, one may 
ignore it. However, as A decreases and the size of the universe grows, radiation becomes 
important. More radiation also means stronger suppression of the tunneling, consistent 
with the intuitive picture of tunneling while interacting with the environment. It is this 
suppression that provides a bound to the Euclidean gravity action. 

7. Connection to SOUP 

In [3] we lay out the motivation for SOUP, short for "Selection of the Original Universe 
Principle". This is an alternate to the Anthropic Principle. Since observational evidence 
of an inflationary epoch is very strong, we suggest that the selection of our particular 
vacuum state follows from the evolution of the inflationary epoch. That is, our particular 
vacuum site in the cosmic landscape must be at the end of a road that an inflationary 
universe will naturally follow. Any vacuum state that cannot be reached by (or connected 
to) an inflationary stage can be ignored in the search of candidate vacua. That is, the 
issue of the selection of our vacuum state becomes the question of the selection of an 
inflationary universe, or the selection of an original universe that eventually evolves to 
an inflationary universe, which then evolves to our universe today. The landscape of 
inflationary states/universes should be much better under control, since the inflationary 
scale is rather close to the string scale. In [3] we proposed that, by analyzing all known 
string vacua and string inflationary scenarios, one may be able to pin down SOUP. Here, 
the rate of tunneling from nothing (i.e., no classical spacetime) to a deSitter universe (much 
like the inflationary universe) is now given by T ~ exp(F) = exp(— SE,dc), where 

-F = S E , dC = S Efl + D^-^ R + n ^ (7.1) 

where D is the decoherence term, l s is the cut-off scale and n^of is the number of light 
degrees of freedom. Here V4 is the 4-volume of the instanton. For S* 4 , V4 = 8-7r 2 r 4 /3 = 
24tt 2 /A 2 is its area. In string theory, that cut-off is naturally provided by the string scale. 
Note that SE,dc is now bounded from below. 
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To gain some idea of the magnitude of F, we use Eq.(|l.S|) to crudely estimate the value 
of A with maximum tunneling probability, 



4n dof G 

A ~ - 

7U„ 



so that 



no 



" (2vrM„/ s ) 4 



2rado/ 

For the string scale a few orders of magnitude smaller than Mpi, we see that A max takes a 
value quite close to that expected in an inflationary universe, with F ~ 10 12 . On the other 
hand, the critical value of A is A c ~ A max /2. For A < A c , the barrel-shaped instanton 
is destroyed. At A c , F(A C ) ~ 0, so the tunneling probability at A close to A c is already 
negligibly small compared to that at A max . Tunneling to a supersymmetric vacuum is 
totally suppressed. For l s around the TeV scale [28], we find that A max is quite close to 
today's dark energy value. This is similar in spirit to Ref [29]. However, this scenario will 
imply that our universe has not really gone through the whole hot big bang and certainly 
not the standard inflationary epoch. 

The key tool is the modified bounce. Although we only deal with the S* 4 instanton and 
its modification in this paper, one can easily generalize this analysis to higher dimensions. 
We expect the form of the tunneling probability to be P ~ e F , with F = —SEdC = 
-S E - D. In 10-D, 

S E ,dC - Se,io + c (7.2) 

where Se,io is the 10-D Euclidean action determined in mini-superspace and Vio is the 
10-dimensional volume of the instanton. In effective 4-D theory, Se io reduces to —3ir/GA, 
since 

8ttG = M- 2 - gs s 



Pl ~ AttV 6 

where Vq is the 6-D compactification volume. (Recall a' = M~ 2 = (I s /2tt) 2 .) The constant 
c has to be calculated for each vacuum and will depend on the details of the vacuum. To 
get an order of magnitude estimate of c we have, by comparing to the 4-D case, 

c ~ nd °f 1 (7 3) 

C vr M 2 pi l 2 g 2 {L6) 

so we expect c to be small. In [3] we have considered such instantons. We did not do a 
careful calculation of v then, that has been the main content of the present paper, but we 
phenomenologically guesstimated a possible range for its values and calculated the creation 
probabilities for various vacua in the cosmic landscape. Among all the known string vacua, 
we find that the universe most likely to be created is a KKLMMT type inflationary vacuum. 
(For details, we refer the reader to [3].) The probabilities calculated were very robust and 
did not depend on the exact value of c. Changing from Vg to V\q does not change the 
overall qualitative results. Therefore, we expect the conclusions to stay the same. 
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There are some minor differences. There we find that tunneling to a universe with 
today's dark energy is very much suppressed. Here we see that the cosmological constant 
corresponding to today's dark energy is actually below the critical value, so tunneling di- 
rectly to today's universe is simply zero. A more careful calculation based on Eq.([T.2j) 
is clearly needed. However, the important point is that one can calculate it, at least in 
principle, for any given vacuum and, therefore, compare their probabilities. This would be 
a huge improvement over the Anthropic Principle. 



One can calculate the tunnelling probability to S 10 using Eq^T^). For S 10 , Se, 10 = 
3^^, where V 10 = ^^/H 10 , is the volume of S w . Note that for S w , H 2 = A/36. 
Without branes, the solution we consider is a 10D super symmetric vacuum. Tunneling 
to this vacuum is suppressed. Next we consider tunneling to a deSitter-like vacuum. Let 
there be N D9 brane-antibrane pairs in the S 10 . So 

M 10 a M 2 N 

Maximizing F in terms of N gives 

N= 5n dof (2,f W» 

9s \MpiJ y ' 

This means A = 1 is a reasonable choice for g s ~ 1, and M s /Mpi = 10~ 3 . One can then 
calculate F for N = 1 for 5 10 . In this case, F is dominated by Se,o- One gets 

F(S W ) ~ -S Efl (S w ) = 5 - g 6?r)5 ~ 4 x 10 9 (7.6) 

Following Ref [3], we see that other related geometries such as S 5 x S 5 , S A x S 6 , 
S 4 x S 3 x S 3 etc. have slightly smaller but similar values of F. 

For a KKLMMT-like inflationary universe, the choice of fluxes fixes the string scale 
as well as the scale of inflation, where G and the density perturbation measured in the 
cosmic microwave background radiation are used as input parameters. Using Eq.([0|) and 
A = 8irGp vac , where the compactifled volume has been absorbed into G, one can rewrite 
Eq.Q so 

3 / V8^M Pl \ 4 3 / M 8 \ 4 ( V8^M Pl V 

Plugging in the above mentioned values for the mass scales, we get F. For a typical choice 
of fluxes that would give M s ~ 10 15 GeV and an inflationary scale plit ~ 10 14 GeV, with 
the Planck mass Mpi ~ 10 18 GeV, one gets 

F ~ 10 16 

for a realistic brane inflationary scenario. This big increase in F from the 10D scenario 
is a result of the decrease in the effective cosmological constant due to the warping of the 
geometry. By varying the RR and the NS-NS fluxes, one can maximize F in a way which 
is not possible in 10D. Of course, this value is sensitive to the details of the model. A more 
careful calculation will be important. 
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8. Discussion 



Let us make some comments here. 

• The loop correction to the Euclidean action has been calculated before [10, 11]. It 
has the form log(A//x), so Se becomes 

Se = ~GA ~* ~GA (1 + ^ GA1 °g( A /^)) 

where (3 measures the number of fields involved. Ref [32] considers the metric per- 
turbation also; however, they then relate it to the one-loop contribution, yielding the 
above result. Note that this loop correction does not solve the boundedness problem 
in the HH wavefunction. Its implications on probability of inflation has also been 
discussed (see [33] and the references therein). 

• Decoherence effects are to be distinguished from the particle creation effects that 
have been discussed in the literature [9]. Particle creation effects arise at the one- 
loop level, while the decoherence effect discussed here is a back-reaction improved 
quantum correction. 

• Back reaction is crucial in obtaining the decoherence term. This reminds one of the 
situation in quantum field theory, where a simple 1-loop correction to the coupling 
takes on new significance when it is renormalization group improved. 

• Integrating out the environment provides a new interaction term for the cosmic scale 
factor a in the Einstein equation. This decoherence term behaves just like ordinary 
radiation and generates an effective term in the Euclidean action for the instanton. 
The presence of this term provides the lower bound to the resulting Euclidean gravity. 

• The correction to the Hartle-Hawking wavefunction results in a decrease in the Eu- 
clidean action. If we are to interpret the Euclidean action as the entropy, then there 
is a decrease in the entropy due to decoherence. This is in accordance with the 
Bckenstein bound on entropy. A pure deSitter space provides the upper bound on 
entropy. Decoherence, which is due to the inclusion of extra degrees of freedom, then 
leads to a lower entropy respecting the entropy bound. 

• Once the inflationary universe is created, inflation simply red-shifts the radiation so 
the radiation term quickly becomes negligible. It seems that the initial radiation is 
present just to suppress the tunneling. 

• It is interesting to speculate what happens after inflation. Towards the end of infla- 
tion, the inflaton field (or the associated tachyon mode) rolls down to the bottom of 
the inflaton potential. The universe is expected to be heated up to start the hot big 
bang epoch. It is not unreasonable to expect the wavefunction of the universe to be 
a linear superposition of many vacua (say, a billion of them) at the foothill of the 
inflaton potential. This is like a Bloch wavefunction as proposed in Ref [35]. Since 
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tunneling happens only between vacua with (semi-)positive cosmological constants, 
this wavefunction should be a superposition of vacua with only positive cosmological 
constants. The ground state energy of this wavefunction is expected to be much 
smaller than the average vacuum energy, a property of Bloch wavefunctions. As the 
universe evolves, tunneling among the vacua will become suppressed, due to cooling 
as well as decoherence. Eventually, the environment collapses the wavefunction to 
a single vacuum. Since the Bloch wavefunction is able to sample many vacua, it is 
natural to expect that it will collapse to the vacuum state with the smallest positive 
vacuum energy within teh sample. This may partly explain why the dark energy is 
so small. 

• In terms of recent work on the statistics of the landscape [30,31], one may view our 
work as providing a measure to the counting of vacua. Each vacuum is weighted 
with its probability of tunneling from nothing. Since the exponent of the tunneling 
probability can differ by many orders of magnitude, this is likely to be the dominant 
contribution to the measure. We argue that a 4-dimensional inflationary universe 
very much like ours has large probability; this implies that vacua that cannot be 
reached after inflation should have vanishing measure. 

• The effect of decoherence on other tunneling problems in quantum gravity should 
be re-examined. These include tunneling in eternal inflation, the Coleman-deLuccia 
tunneling etc. 

9. Summary and Remarks 

We address a number of questions in this paper. 

First we ask what happens to the well known S 4 bounce solution to the Euclidean 
Einstein equations when the perturbations to the metric are taken into account. How 
does the bounce change as a result? In particular, is there a parameter characterizing 
the perturbations that describes the modification of the bounce? Is there a range of 
this parameter that destroys the bounce altogether? 

To answer the first question, we apply the path integral techniques to trace out the 
perturbations. The result is summarized in the introduction. Taking back-reaction 
into account, the Euclidean action now includes an additional term that encapsulates 
the effect of the perturbations. Depending on the value of the parameter c or equiv- 
alently v, the bounce gets deformed. For a fairly large range of this parameter, the 
effect is to just deform the bounce solution. Since this contribution is always positive, 
the tunneling is always suppressed. At a critical value, the bounce is destroyed. One 
may interpret that tunneling is forbidden in this limit. 

Next, we ask what these perturbations do to the boundedness of the Euclidean grav- 
itational action. As is well known, the Euclidean gravitational action for a closed 
spacetime is not bounded from below for theories with a dynamical cosmological 
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constant. This manifests, for example, as the infinite peaking of the wavefunction 
of the universe at the vanishing value of A. Does the inclusion of the perturbations 
change anything here? 

To answer this question, we see that the effect of the perturbations is to change 
the wavefunction of the universe from exp(3vr/GA) to exp(37r/GA - C/A 2 ) where c 
depends on the particular perturbations that we are looking at. So, at least for the 
problem at hand, the inclusion of the perturbations and their back reaction makes the 
Euclidean action bounded from below. This renders the wavefunction normalizable. 
In the study of tunneling, it makes no sense to talk about gravity without taking into 
account the perturbations to the metric. In this sense, quantum gravity is consistent 
as long as we are careful to include the effects of the perturbative modes which are 
always present. 

Once we have a sensible wavefunction, we can now go ahead and apply it to the 
cosmic landscape in string theory. Since Euclidean action and tunneling probability 
are both dimensionless, we can compare the tunneling from nothing to any point 
in the landscape and find the sites that have the largest probability. This program 
hopefully will allow us to understand why we are where we are, without resorting to 
the anthropic principle. 
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A. A Quantum Mechanical Example 

To see the basic idea, recall the quantum tunneling of the system 

Lo = ^q 2 - V(q) (A.l) 

with a quartic V(q) as shown in Fig. 2. The tunneling rate T from the local minimum 
at q = to the exit point q = qo is well-known, 

r = Aexp(-So) 

S = J" ^2MV(q)dq = J (^-q 2 + V(q)j dr (A.2) 
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where r is the Euclidean time and Sq is the bounce, i.e., the instanton solution [13]. 
This WKB approximation is good provided that the height of the barrier is larger 
than wo, where 

MuJn 



> d 2 v, 



Note that, for V(q) bounded from below, So is bounded from below, as required by 
consistency. In a more realistic situation, the particle interacts with the environment. 
Typically, this introduces a frictional force, so the corresponding classical equation is 
given by 

dV 

Mq + rjq+ — = (A.3) 

The impact of such a frictional term on the particle is to suppress the tunneling rate. 
Consider the following system 

L = y q 2 ~ V(q) + ~ £ m a {x 2 a - J&) - q £ C a x a - ~M(6u) 2 q 2 (A.4) 

a a 

where one may consider q to be the system and the x a to be the environment. The 
last term is a counter term introduced to correct the shift in frequency, 

The interactions of q with the x a introduces the friction term 

for uj smaller than some critical uj c . 

The tunneling rate of this sytem can be easily found [14, 15] that the bounce So 
increases to 

s = s »( i+ 2]>y < a - 7 > 

That is, the interaction with the environment, or the friction, suppresses the tunneling 
rate. This qualitative feature remains true when the environment and its interaction 
with the system is more complicated. 



B. Calculations on the Squashed S* 4 

This appendix gives details for Section 3. We show how we perform the path integral 
over the t n mode. The solution to Eq.(|3.1|) satisfying the initial condition t n (x = 
1 — 5) = t l n , and the final condition t n (x = —1 + 5) = th, is given by 

■ f , _ F n (x) AP 1 - n (x)+BQUx) . . 
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where A and B are given by 
A = y/25 - 5 2 — 



t j n Ql{-l + S)-tlQT[(l-5) 



B = -^25 - 5 2 — 



Pf"(l - 5)QI{-1 + 5)- Pf n (-1 + 5)QU1 - S) 

tip- n (-i+d)-tipr(i-s) 



(B.2) 



Pf n (l - (-1 + 5) - Pf n (-l + " -5) 

The trace operation is defined as doing the following 

(j-initial xi \ 

V L n ~ L n) 



< / dtfj{t - tl) 



fj.final_.fs 
\ L n — L n ) 



D[t n ]exp (£%[t n (x)]) 



(B.3) 



Since our purpose finally is to take to trace over the t n mode, we can set t l n = tn at 
this stage. The Euclidean action due to the nth mode is then 



5n 
w. 



a{rft r 



dt n 
l dr 



x=(-l+S) 
x=(l-5) 



(B.4) 



xF 2 + (1 - x 2 )F r 



dF n 
dx 



x=(-l+S) 



x=(l-S) 



Using Eq. i 



B.2 ) and the symmetry properties of the associate Legendre functions 
Pr n (~x) = -Pr n (x), QU-x) = QUx) (B.5) 



d x P{ n (-x) = d x P{ n (-x), d x Q n x (-x) = -d x QUx) 
one gets the following action 



S n E = 5(2 - 5) 



(1-5) + 5(2 -5) 



d x Q1(l - 5) 



(t 



i \2 



(B.6) 



Q?(l - S) 

Next we must find the prefactor to the path integral for t n . This can be found as 

where g n (x) 

CPi n (x)+DQ1(x) 



explained in the Appendix C. The prefactor is given by , 1 = 



is a solution to Eq, Q3.1| ), i-e. g n (x 
conditions 



VT- 



and it satisfies the following 



g n (l-5) = 0, d x9n (l-5) = l (B.7) 

It is easy to check that such a solution to Eq.([0|) has the following values of C and 
D 



C= ^25 - 5 2 
D = -V25 - 5 2 - 



Q?(l - 5)d x P{ n (l -5)- Pf n (l - 5)d x Ql(l - 5) 

P__ ~ S) 

QUI - 5)d x P{ n (l -5)- Pf n (l - S)d x Cft(l - 5) 
The prefactor is then given by 

1 1 



(B.8) 



8 x Q^(l-5) d x P{ n (l-5) 

^27rg n (-l + 5) v^'V V " $) ^f n (l - 5) 



(B.9) 



i v W r WMM /r(3^)r(2 5 a) 



4tt 



2" 
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C. Path Integral of an Oscillator with variable mass and variable 
frequency 



Here we review some basic properties of path integral and apply them to the evalua- 
tion of the prefactor in Eq.(4.7). To be concrete, we shall follow Ref [34]. Consider an 



oscillator with mass m(t) and frequency u(t). What we have in mind, in particular, is 



a case like Eq.(4.2) which describes the path integral for a tensor mode (in Euclidean 
spacetime). Comparing Eq.([4.2|) with this appendix, the time dependent mass would 
be a(t) 3 and the time dependent frequency would be (n 2 — l)/a 2 . We would like to 
evaluate the propagator K (xf, T\xi, 0). The action is given by 

S = £dt(\m{^) 2 -\muj(t?Xz) (C.l) 

The problem can be simplified by introducing a new "time" variable u so as to map 
the present problem to that of an oscillator with unit mass and variable frequency: 

du = -4t (C.2) 
m(t) 

In terms of u the action becomes 

s = \Qdu([^) 2 -n{ufx^) (c.3) 

that is, the action for an oscillator with unit mass and variable frequency Q(u). As 
this is quadratic in X(u), we can expand around a classical solution, X = x c i + x, 
where the classical solution x c i (u) satisfies the equation of motion 

^^ + {l{ufx d {u) = V (C.4) 

where Q(u) depends on n. For example, for the nth tensor field perturbation, Q(u) = 
(n 2 — l)a('u) 4 . One can use the above solution for x c i, where x c i{ui) = X(m) and 
x d( u f) = X(uj), to calculate the classical action S^a;^]. This gives the saddle point 
value of the path integral. Following the discussion in Sec. 5, the prefactor of the 
propagator is given by: 

r&5=? D \?m ^ [! c du ((© 2 - ^) 2 * 2 )] ( c - 5 ) 

By a further change of variables, we can change the action to a free particle action. 
To do so, let g(u) be a solution of the equation of motion (|C.4| ) such that g(u) does 
not vanish at the initial end point 

g{ Ui ) ± (C.6) 

implying that g(u) is not a path in Eq.( |C.5| ). Define the following transformation of 
variables 

x(u) = g(u) P ^ds = g(u)F(u) (C.7) 
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where the function y(u) obeys y(ui) = 0. Differentiating Eq.( |C.7| ), one finds the 
inverse transformation 



d 2 
dii 2 



ds 



y(u) = x(u) - Q ^-x(s)ds 
In terms of the y variable, we have: 

+ n(u) 2 ) x(u) = (g"(u) + n(u) 2 g(u)) £ 

I g'("tf(") i -j/'tnA — 9'{u)y'{u) /// 

so the action becomes 

S[x(u)] = -\ HI du [F(u)g'(u)y'(u) + F(u)g(u)y"{u)] 



9(s) 



(C.8) 



(C.9) 



(CIO) 



Keeping in mind that x(u) vanishes at the boundaries (see Eq.( |C.5| )), a further partial 
integration then leads to 



This is just a free particle action, with boundary conditions 



y{ui) = ; y(uf) 



Lf ds m=o 



g(s) 



(C.ll) 



(C.12) 



To impose the (non-local) boundary condition x(uf) = y{uf) = 0, we introduce 

5(x(u f )) = ±fda e~ iax( - u f^ (C.13) 
The path integral can now be written as 

D[x{u)] exp { iS[x(u)]) = £ CSl^ fZ daD[y(u)] 



det 



5x 




Sy 


. exp 



-iag(u f ) Q ds^- { 



. exp 



2 -in; aa \ du 



(C.14) 



As the transformation between x(u) and y(u) is linear, the Jacobian det 
independent of y(u). To carry out the above path integral, we just have to comp 
the square by the use of the new variable 



is 
ete 



7(u) = y(u) - aug(u f ) Q ^ 



The above path integral thus becomes 



2tt 



det 



fy(uf)arbitrary 



JZ da exp y-L^gtuf) 2 Q 

2 



(C.15) 



(C.16) 



D[y(u)] exp 



1 r u f ( gh\ 

2 Jui \du J 
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It is easy to carry out the a integral. To do the remaining path integral, one just 
has to notice that it represents the probability amplitude for finding the particle 
anywhere at the time Uf. This probability is obviously 1. 



= J^° oo dxK(x,Uf\0,u i ) = 1 
where 7(w/) is arbitrary. Therefore, we have 



r-x('L 

Jxlu 



det 



(C.17) 



(C.18) 



It remains to find the value of the Jacobian. It can be found by a discretization 
process. The paths x(u) and g{u) can be replaced by the multidimensional points 
(x ,x 1 ,...,xn) and (y , y%, y N ) with x k = x(u k ) and y k = y(u k ). The linear 
transformation can then be approximated as 



Then, 



Vn 



X 



_ T g'(u k ) x k + x fc -i 

n n hi ^ uk) 2 



(C.19) 



Jn = det 



% 
dxj 



N 



n 1-5 



1 5'(«fc) 



fe=l 



Taking the N — > oo limit, one gets 



det 



5y_ 

8x 



exp 



So the path integral prefactor becomes 

rx(uf)=0 



A 



1/2 



S(u f ) 



n( u ) 2 x 2 



det 



du 2 



n( u y 



-1/2 



(C.20) 



(C.21) 



(C.22) 



where A is a normalization factor. One can check the validity of this result for 
simple cases like free particle and simple harmonic oscillator. Using the result from 
Eq.( |C.22"D , we can write 



det 



. d 2 
du' 2 



g( Ui )g(u f )£ 



u f ds 



det 



d 2 

du 2 



u{u) )] G( Ui )G( Uf ) /;/ g*. 



(C.23) 



- 34 - 



where G(u), analogous to g(u), is a solution of the equation of motion with an 
arbitrary U(u). And, like g(u), G(u) also does not vanish at Uj. Now, let g° denote 
the unique solution to + Q(u) 2 g(u) = which satisfies the boundary conditions: 



g°(m) = 0; = i 



Let g 1 represent the solution that satisfies 



In Eq.(^]2|), let 



9 1 (u l ) = l; —g 1 (u i ) = 



g(u) = g°(u) + eg^u); G(u) = G°(u) + eG 1 ^) 



(C.24) 



(C.25) 



(C.26) 



The integral f™ f \J%± diverges due to the vanishing of g° at U{. However, since 
almost all the contribution comes from an infinitely small neighborhood of Ui (in the 
limit when e — ► 0), it follows that it diverges like 



du 



U - Ui 



Consequently, 



lim 

Ju 



So Eq. (|C.23|) becomes 



det 



. d 2 
du 1 



Uf du 



n(«) 



det 



<P 



du* U ( u 

Taking U(u) = (for a free particle), we get 

det (-£-n(u)* 



ruj du 
Jui Giu) 1 



9>f) 
G°(u f ) 



(C.27) 



(C.28) 



(C.29) 



det 



d? 

du 2 



9>f) 

Uf - Ui 



(C.30) 



Now we can calculate the propagator, 



K(xf,T\xi,0) = K(xf,Uf\xi,Ui) = K(0,Uf\0,Ui) exp (iS[xci\) 

-1/2 

- K (0,Uf\0,Ui)exp(iS[x d ]) 



dot —2^ -Q(uy 



det - 



(C.31) 



One can switch back to the t variable now using Eq, flC.2|) . 5[a; c /] is the action for 
the solution x c i of the classical equation of motion satisfying the boundary conditions 
x(ui) = Xi and x(uf) = Xf. 
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Now we specialize this calculation of the path integral for an oscillator with varying 
mass and frequency to our case. Eq.( p,4j) generally does not lend itself to simple 
analysis. However, since we expect the higher modes (with large values of n) to be 
the major contributors to the suppression of the tunnelling, one can use the WKB 
analysis to find approximate solutions. This is because, firstly, Eq.( |C.4| ) has the 
form of a Schrodinger equation for an arbitrary potential Q(u) 2 , and, secondly, the 
variation in £l(u) 2 x c i(u) for higher modes is slow enough. Typically, for the nth mode, 



fi 2 



n 



2 and Cl ~ n. So Cl « Q 2 . Using the WKB approximation, we get: 



Xd(u) 



exp ±i 



duCl(u) 



(C.32) 



Note that in the Euclidean case, u — > —iu (or, equivalently, t — * —it). This is what we 
have done in Eq.( f4.1C , 4.11 ) to get the solution for the tensor mode t c i(r). Using these 
classical solutions we have calculated the classical action in Eq.( [4.15| ). As explained 
in Eq.(C.31), the prefactor is proportional to {2 / nig a {uf))~ 1 / 2 where g°(u) satisfies 
the boundary conditions given in Eq. (|C.24|) . Following Eq.(gl|) it is clear that 
for tensor modes, the two linearly independent solutions to the classical equation of 
motion are exp (± J u du'uj n (u')^ . And the following linear combination satisfies 
the required boundary conditions 



9°(u f ) 



yj2io n (ui)u n {uf) 



(exp(D„) - exp (-£>„)) . 



(C.33) 



as this choice of g (it) vanishes at u 
yields the result in Eq.( [4.16 ). 



Ui where D n = 0. 



dg°{u) 
du 



1 at « 



Ui. 



This 
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